Collective modes of trapped gases at the BEC-BCS crossover 
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The collective mode frequencies in isotropic and deformed traps are calculated for general polytropic 
equation of states, P oc n^^^, and expressed in terms of 7 and the trap geometry. For molecular and 
standard Bose-Einstein condensates and Fermi geises near Feshbach resonances, the effective power 
7 ~ 0.5 — 1.3 is calculated from Jastrow type wave-function ansatze, and from the crossover model 
of Leggett. The resulting mode frequencies are calculated for these phases around the BCS-BEC 
crossover. 



Recent experiments probe systems of fermions [1-5] 
and bosons [6] near Feshbach resonances by expansion 
and RF spectroscopy. Interesting new strongly interact- 
ing or dense phases of bosons and fermions are created, 
e.g., that associated with the crossover from a (possibly 
superfluid) Fermi gas to a molecular BEC. The corre- 
sponding equations of states (EOS) differ from standard 
dilute systems which directly shows up in their collective 
modes. 

The purpose of this work is to calculate the collec- 
tive modes in terms of a general class of polytropic EOS, 
to calculate the EOS for strongly interacting BEC and 
Fermi gas at the crossover to a molecular BEC, and fi- 
nally synthesize the two to calculate the collective modes 
for these strongly interacting phases. 

Polytropic EOS relate the pressure and density as 



P cx n''+^ . 



(1) 



As we shall show below the collective modes in harmonic 
oscillator traps depend on the power 7 but not on other 
details of the polytropic EOS. Polytropic EOS apply to 
many systems. In a dilute interaction dominated BEC 
7 = 1, whereas an ideal Bose gas in the normal state has 
7 = 2/3 under adiabatic conditions. A dilute gas of Fermi 
atoms also has 7 = 2/3 in both the hydrodynamic and 
superfluid limits. Both a Fermi gas [7] and a BEC [8] has 
7 = 2/3 in the strongly interacting (unitarity) limit. We 
shall see below that near a Feshbach resonance, where 
a Fermi gas crossover to a molecular BEC, the power 
effectively varies between 7 ~ 0.5 — 1.3. 

The collective modes are calculated from the equations 
of motion which in hydrodynamics and for a superfluid 
are given by the equation of continuity and the Euler 
equation 



mn— = -VP - nWext ■ 
ot 



(2) 



Here n and v are the local density and velocity, and 
Vext = (l/2)m^j w^rf is the harmonic oscillator trap 
potential. From the Euler equation we obtain the equi- 



librium density: 



"0(1 - T,,r!/RlV^'', where 
i?2 = 2{j+l)Po/jnomu;f, i = 1, 2, 3, are the 3D Thomas- 
Fermi radii of the trapped cloud of atoms (Pq and uq are 
the pressure and density in the center of the trap) . 



Linearizing around equilibrium, n = n^„ + e^'^^Sn, the 



equations of motion lead to 
—muP'Sn — V 



nV I —dn 

n dn 



(3) 



It is not necessary to restrict ourselves to zero tempera- 
ture where the Gibbs-Duhem relation dP — ndfi simpli- 
fies Eq. (3). 

In an isotropic trap the collective modes with angu- 
lar momentum I and n radial nodes are straight forward 
to calculate from Eq. (3) by generalizing the method of 
Ref. [9] to any polytropic EOS. We find that the de- 
parture from the equilibrium density is Sn oc r'(l — 
rVp2)(i/7-i)_p(_„^„ + I + 7-1,/ + 3/2,r^/R^), where 
F is the hypergeometrical function. The corresponding 
eigenvalues are 



=; + 2n[7(n + / + l/2) + l] 



(4) 



which reduces to the known results for 7 = 1 [10] and 
7 = 2/3 [11]. In comparison the collective modes in the 
collisionless limit arc those of a free particle: uj/ujo — 
2n + I, when its mean free path exceeds the size of the 
cloud. 

The hydrodynamic collective modes can also be calcu- 
lated for deformed traps for a general polytropic EOS. 
Linearizing the equations of motion lead to the following 
equation for the collective modes 



-LU^V = V(V • VVext) + jC^Vext)C^ ' v) . 



(5) 



The breathing modes have flow velocity on the form 
V = (ttiri, a27'2, a3r3)e"^*, which leads to three coupled 
homogeneous equation for ai=i,2,3 



^Y.'^ho = (w' - (7 + 2)^,?; 



(6) 



In an axial symmetric trap: lui — lu2 = luq and lu^ — Xluq, 
the resulting breathing modes are [12] 



^ =7 + 1 



7 + 2 



± V(7 + 2)2AV4 + (72 - 37 - 2)A2 -I- (7 -I- 1)2 , (7) 
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and cij = \/2^o- The ± eigenvalues are the radial and 
axial modes respectively and result from the coupled 
monopole and quadrupole m = modes, where m is the 
angular momentum projection on the 3rd axis. There- 
fore, the breathing modes for an isotropic trap A = 1 be- 
come the quadrupole with oj ~ V2uJo and the monopole 
with ui = y'37 -|- 2cl)o as follows from both Eq. (7) and 
Eq. (4) for n = 0,Z = 2 and n = 1,Z = respectively. 
For 7 = 1 and 7 = 2/3 the monopole frequencies are the 
standard u> = V^loq and to = 2wo respectively. 

For a very elongated or cigar-shaped trap (prolate in 
nuclear terminology), A <C 1, used in recent experiments 
[5], Eq. (7) results in a low frequency axial mode with 



= V3-(7 + l)- 



U)5 . 



(8) 



For 7=1 and 7 = 2/3 the axial mode frequencies are 
w = \/5/2uj3 and u) 



^12/5a;3 respectively as found in 
[13,14]. The radial modes have 



(^rad 



^2(7+1) WO . (9) 

For 7=1 and 7 = 2/3 the radial mode frequencies are 
the standard uj = 2uja and lo = ■y/lO/3a;o respectively. 

In the oblate limit, A ^ 1, the breathing modes are 
Wax = V7 + 2a;3 and uJrad = V(67 + 4)/(7 + 2)wo. They 
connect to the prolate limit through avoided level cross- 
ing as seen in Fig. 1. 
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FIG. 1. The ± solutions of Eq. (7) vs. trap deformation 
A = ujs/loo for 7 = 1/3, 2/3, 1, 4/3. For a cigar shaped trap, 
A <C 1, these correspond to the radial and axial mode fre- 
quencies respectively. 

We now turn to the EOS for strongly interacting Bose 
and Fermi gases and calculate an effective polytropic 
index that can be applied for the above modes. The 
EOS for a BEG was calculated in Ref. [8] from a Jas- 
trow type wave function ^'./(ri, ...,rAr) = Y[i<i 
which incorporates essential two-body correlations and 
is a good approximation for cold dilute and dense bose 
systems [15]. It was shown that with proper bound- 
ary conditions the calculated energy reproduced the di- 
lute limit result, E/N = 2'Kh'^an/m, where a is the 



s-wavc scattering length between bosons. However, in 
the imitarity limit, n^^^a ^ 1 or x = 1/akp ~ 0, the 
energy per particle scales like a Fermi gas polytrope: 
E/N = 13.33n'^n'^/^/m = 2.79Ef. Here, we have also 
for bosons defined Ep = T)?k%/2m in terms of the den- 
sity n = kp/Sn'^ for later comparison between molccnilar 
BEG and Fermi gases with two spin states. In this model 
(see [8,16] and Fig. 2 for details) we can calculate the zero 
temperature pressure P = d{E / N) j dn and the effec- 
tive polytropic index, which we define as the logarithmic 
derivative 



7 = 



n dP 



- 1. 



(10) 



As shown in Fig. 3 7 approaches 1 asymptotically from 
above in the dilute limit but is 2/3 in the unitarity limit 
with a maximum of ~1.25 around l/kpam ~ 1-4. 
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FIG. 2. The energy per particle in units of Ep for a BEG 
with ttm = a and am = 0.6a, and a Fermi gas as it crossover 
towards a molecular BEG (see text). At small positive scatter- 
ing length the energy per fermion approaches half the binding 
energy of a molecule, E/N —h^/2ma^ = —x^Ep. 

The approximate EOS can also be applied to a molecu- 
lar BEG when a number of factors are taken into account: 
the density of molecules is half that of Fermi atoms, their 
mass is two times larger, the above calculated E/N is 
for two atoms, and finally the scattering length between 
molecules (a^) may differ from that between the two 
Fermi atoms forming the molecule. Petrov et al. find 
[17] Qm — 0.6a in accordance with [3]. As a result the 
energy in molecular BEG is per atom much smaller than 
that of a BEG of Bose atoms as seen in Fig. 2. This 
is important when we attempt to match the EOS of a 
molecular BEG onto a gas of Fermi atoms in the unitar- 
ity limit. The difference in magnitude cancels in Eq. (10) 
and is therefore not important for calculating 7. A dif- 
ference between a and and between Bose and Fermi 
densities does, however, affects akp as seen in Fig. 3. 

The EOS for a Fermi gas at zero temperature near 
a Feshbach resonance and its crossover to a molecular 
BEG has recently been studied by various resummation 
techniques [18,7], a Jastrow-Slater type ansatz [7,16] and 
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by fixed-node Greens function Monte Carlo (FN-GFMC) 
[19]. The EOS calculated from the Jastrow-Slater ansatz 
has the merit that the EOS is exact to leading orders 
both in the dilute and molecular limits. Furthermore, it 
has been tested experimentally [16] and in FN-GFMC 
to be a good approximation in the unitarity limit as 
well. It extends the Jastrow wave function for bosons 
described above by including an anti-symmetric Slater 
wave function ($s), which is the ground state of free 
fermions *js(ri, tat) = '^sU.i^j^ fi^i-^j')- Because 
$5 insures that same spins are spatially anti-symmetric, 
the Jastrow wave function only applies to particles with 
different spins (indicated by the primes). The resulting 
energy (see [19,16] for details) is shown in Fig. 2. The 
energy per particle at zero temperature can generally be 
written in terms of the ratio between the interaction and 
kinetic energies /? = Eint/Ekin as [7,1] 

E/N = Ekin + Eint = ^Ef[1 + 13] . (11) 

As seen in Fig. 2 it approaches the ideal Fermi gas result 
(3/5)^;^ in the dilute limit, (3/5)(l + (3)Ef in the uni- 
tarity limit with (3 ~ —0.54 for two spin states [16], and 
E/N = ~h^/2ma?- in the molecular BEC limit (a ^ 0+). 
The corresponding pressure is 

P = n'^^J^ = lEFn[l+l3-xl3'/2]. (12) 
an 5 

Here, we view /?(.x) as a smooth function of x = \/akp 
with derivatives (3' = d(3{x)/dx, etc. The effective poly- 
tropic index is from Eqs. (10) and (12) 

__ |(l+/?)-x/?72 + xW6 

^- l + l3-xl3'/2 • ^^"^^ 

When I3{x) is a smooth function at x = we find from 
Eq. (13) that 7 = 2/3 in the unitarity limit. 

The effective polytropic index is shown in Fig. 3 for 
a Fermi gas as it crossover to a molecular BEC. That it 
turns over and drops back to 7 ^ 2/3 for a;^0.5 is an 
artifact of the EOS resulting from the Slater ansatz in 
the wave function. The true ground state wave function 
is expected to have a lower energy in the molecular BEC 
limit as a ^ as is also found in FN-GFMC [19]. The 
Jastrow part of the wave function is responsible for the 
correct leading part of the energy: E/N = —h^/2ma? 
which, however, does not contribute to the pressure be- 
cause it is density independent. The Slater wave func- 
tion is responsible for the leading density dependent or- 
der but detailed comparison to FN-GFMC calculations 
show that it is only correct up to a;;S0.5. In FN-GFMC 
a better ground state wave function is found numerically 
which has lower energy. Both Fermi gases and BEC's 
have 7 = 2/3 at x = due to the universal scaling law 
E/N oc n^/^ in the unitarity limit [7,8]. 

Leggctt [20] extended the BCS gap equation to the 
BCS-BEC crossover and calculated the gap and chemi- 
cal potential. Prom the Gibbs-Duhem relation and Eq. 



(10) we can then calculate 7 as shown in Fig. 3. In the 
dilute BCS limit it differs from the Fermi gas, which has 
chemical potential n — Ep + 2'Kfi' an/m, by lacking this 
second term proportional to a. In the other limit the 
chemical potentials differ for orders higher than linear in 
the scattering length. In both the dilute BEC and the 
Jastrow BEC approximation higher orders add positively 
- which is responsible for 7 > 1 for fcirOm^l. The Legget 
model leads to negative higher order contributions be- 
yond the linear one: yit = — ?i^/2ma^ -|- ■nh^an/m, and 
therefore 7 decreases monotonously from 1 towards 2/3 
as a; ^ 0+. 

In all the above models the two-body correlation func- 
tion undergoes a smooth transition from a constant in the 
dilute Fermi gas limit to that of a molecule in the dilute 
molecular Bose gas limit, and the BCS-BEC crossover is 
continuous. In the unitarity region the correlation length 
is of order the interparticle spacing r-^ kp^; only in the 
dilute BEC limit a; 1 is the correlation length suffi- 
ciently small that the molecules may be approximated as 
point particles. 
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FIG. 3. The polytropic power 7 of a BEC with am = a 
and dm ~ 0.6a, a Fermi gas (Jastrow-Slatcr ansatz), and the 
Leggett model describing the smooth crossover to a molecular 
BEC. See text and Fig. 2. 

The strongly interacting EOS's near the unitarity lim- 
its can be approximated by a polytrope by replacing 7 
with 7. This allows us to calculate the collective modes 
directly from Eqs. (4) and (7) as shown in Fig. 4. Ex- 
perimentally one tunes the scattering length near a Fesh- 
bach resonance for a fixed number of trapped particles A'^ 
whereby the size and density of the cloud and therefore 
also kp varies in a complicated way depending on the 
EOS. For fermions in the dilute and unitarity limit 7 = 
2/3 and the size of cloud is i? = {24:NY^^aosc{l + 
with = and /? ~ —0.56 respectively. Since 7 ~ 2/3 
for a Fermi gas up to and around the unitarity limit this 
relation for the size is a good approximation in this re- 
gion, and analogously: kp ~ (247V)i/6a-i (1 + /3)V4. 

In a dilute BEC: R = {15 Na„,y^^atic and kp = 
(37r/8)i/3(15iV)2/i5a-i/5a-4/5 in the center of the trap. 
The resulting collective modes - in particular the radial 
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breathing frequency - are as shown in Fig. 4 very sen- 
sitive to a; = 1/akp through 7 near the unitarity limit, 
where the Fermi gas crossover to a molecular BEC. The 
modes can therefore be exploited to extract the EOS ex- 
perimentally. The leading corrections in the dilute Fermi 
gas and BEC were given in [21]. It is observed from Fig. 
4 that the frequencies are very sensitive to am /a, and may 
therefore be exploited to relate the atomic and molecular 
scattering lengths. 
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X = 1/akp 

FIG. 4. The radial and axial frequencies for cigar shaped 
traps, A <C 1, from Eqs. (8) and (9) with 7 from from Fig. 3. 

In very dilute Fermi gases, where the pairing gap be- 
comes smaller than the oscillator frequency and the co- 
herence length exceeds the system size, particle excita- 
tions of order ~ 2 A appear [11,22] besides the collec- 
tive modes described above. Using the pairing gap in an 
isotropic dilute trap [22] the condition A < tiujQ becomes 
X = l/aA;ir<-(2/7r)(C-Mn(3iV)/3), where C = 0.577.. is 
Eulers constant. Such pair excitation modes are therefore 
only observable for weak attractions. 

The collective mode frequencies do not distinguish 
between a superfluid and a hydrodynamic Fermi gas. 
The damping of the modes should be different but has 
not been estimated in the unitarity limit for bosons 
or fermions. Assuming an unitarity limited scattering 
cross section we expect the collision rate to decrease as 
~ exp{—A{T)/T) at temperatures well below the gap 
A(T = 0) ~ 0.54£'i;^exp(7rx/2) [19,16] in a bulk system. 
The damping can potentially discriminate between hy- 
drodynamic and superfluid Fermi gases, and at the same 
time the mode frequency discriminates coUisionlcss. 

In summary, the dependence of collective mode fre- 
quencies and damping on density, interaction strength 
and temperature as described above can - especially near 
the BCS-BEC crossover - reveal the underlying EOS in 
detail including possible phase transitions and associated 
critical temperatures and densities. 

Note added in proof. Two experiments have recently 
published the axial and radial modes around the unitarity 
limit. The above results for the Leggett model is in nice 
agreement with the data of Kinast et al. [23] and (for the 
axial mode) with Bartenstein et al. [24] . See also [25] for 



details. The Leggett model also explains "surprise one" 
in [24] 
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